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ABSTRACT 

We present the first spectroscopic observations in the interval 0.8/im— comple- 
mented with existing HST UV and optical spectroscopy, of the narrow-line Seyfert 1 galaxy 
Mrk 766. The NIR spectrum is characterized by numerous permitted lines of H I, He I, He II 
and Fe II, and forbidden lines of [S ll], [S III] and [Fe ll] among others. High ionized species 
such as [Si IX], [Si x], [S IX] and [Mg Vll] are also observed. The continuum emission has a 
complex shape, with contribution of the central engine, circumnuclear stellar population and 
dust. This last component is evidenced by the presence of an excess of emission peaking at 
2.25/im, well fitted by blackbody function with rbb=1200 K. That temperature is close to 
the evaporation temperature of graphite grains. As such, it provides strong evidence of hot 
dust, probably very close to the nucleus. Consistent modeling of the line spectrum and the 
broad band continuum by composite models, which account for the photoionizing flux of the 
central engine and shocks, shows that shock velocities range between 100 and 500 kms~^, 
the preshock densities between 100 and 1000 cm^^ and the radiation fluxes from the active 
centre between 10^ and 5 10^^ photons cm^^ s~^ eV~^ at 1 Ryd with spectral indices ajjv 
= -1.5 and ax=-0.4. Adopting silicate grains, dust-to-gas ratios are between 10^® and 4 10^^ 
by mass. The emitting clouds are at an average distance of 160 pc from the centre, with high 
velocity clouds closer and low velocity clouds farther from the centre. The N/H relative abun- 
dance deduced from the fit of the [N ll] 6548h-/[0 III] 5007 H- line ratio could be twice solar On 
the other hand, Fe is depleted from the gaseous phase by a factor > 2, most probably trapped 
into grains. Ratios of calculated to observed line ratios to H/3 indicate an average contribution 
of the broad line region to the observed H/3 of about 40%. 
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1 INTRODUCTION 

A large amount of information about the structure and physical con- 
ditions of active galactic nuclei (AGN) has been obtained during the 
last decade from the a nalysis of their opt ical emission-line spec- 
trum (see the review o f|v eron-Cettv & Ve ron 2000 and references 
therein). More recently, spectroscopic observations in the near- and 
mid-infrared have been added to that of the optical with enough res- 
olution to reveal important aspects from both the narrow line and 
broad line regions (NLR and BLR, respectively: [ Sosa-Brito et alJ 
l200lt iKnop et alj|200lt rRodrfguez-Ardila et aljEoOlalbh as wel l 
as from the circumnuclear environment jXhomley et all I2OO0I: 
iLutz et aDl2003j Verma et al. 2003) . Given that the infrared re- 
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gion is substantially less affected by extinction, it also allows 
one to probe depths unreachable with data at shorter wavelengths 
(Sosa-Brito e t al, 2001; Nagar et al. 2002V 

Among the differ ent types of active galaxies, the narrow-line 
Seyfert 1 class (NLSl. lOsterbrock & Pogg? 'l985^ are special tar- 
gets because they share all optical properties of the Seyfert Is ex- 
cept the very narrow Balmer lines (FWHM<2000 km s~^, by def- 
inition) and the strong optical Fe II emission. In the X-rays region, 
NLSl are also peculiar because of the soft X-rays excess above the 
hard X-rays power-law extrapolation that dominates below 1 keV 
I Puchnarewic z et al, 1992 ; Boller_e t al. 199 6) and the very steep 
soft photon spectral index iLeighlviri'999bt) . These extreme prop- 
erties could derive from the fact of NLS 1 being sources with rela- 
tively small black hole masses ( Mbh '"-10^"^ M0) and high E d- 
dington ratios L/Ledd --1-10 (e.g lBorosorj j2002h ; lGruDd j2004b ). 

In the near-infrared (NIR), NL Sl usually present a 
rich emission-line spe ctrum feodrigue z-Ardila et alJ l2002albt 
iGrupe & Thomas! l2002h . In addition, features associated to the 
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BLR such as the Fe II and the permitted H I lines, strongly blended 
in classical Seyfert 1 (Sy 1), are easily separated in NLSl because 
of the narrowness of their BLR lines and the larger separation in 
wavelength between the individual features. This, in turn, favors 
the analysis of the forbidden line spectrum. Weak features, usually 
lost amid strongly blended BLR lines in the optical region, are more 
easily seen in the NIR spectrum. It is then interesting to study the 
cases where the NLR lines are conspicuous in order to examine the 
physical conditions of the emitting gas and the continuum prop- 
erties to construct a coherent picture that explains the parameters 
driving the emission of the NLR in NLS 1 . 

With the above in mind, a program to consistently study a 
selected sample of NLSl galaxies, with emphasis on the NIR in- 
formation, was started with Ark 564 i Contini et al...200 3a: here- 
after Paper I). In this work, a similar analysis is applied to the 
NLSl galaxy Mrk766(=NGC 4253). Located at z = 0.0129, it is 
classified as_a barred SBa galaxy with optical overall orientation 
RA=69° <MacKentv 1990). HST images of this object show fila- 
ments, wisps and irregular dust lanes around an unresolved nucleus 
iMalkan et al. 1998). Radio observation s at 3.6 cm, 6 cm and 20 cm 
lUlvestad et al. 1995VNa gar etalll999t) show that the radio source 
appears to be extended in both P.A. ^27° (on a scale of 0".25) 
and RA. 160° (on a scale of 0"3,'Naaar et al. 1999). In t he op- 
tical, the emission is extended iGonzalez Delgado & Pere'zlll996t 
[Mulchae vetalJll996h through a region of total size greater than 
that of the radio source. 

In X-rays, Mrk 766 is variable on a few hours time- s cale and 
presents a strong soft X - ray excess jMolen di et alj[T993l : iLeighlvl 
ll996tlMason et all2003h . lBoller et alJJl99 ^ report one of the flat- 
test ROSAT soft X-ray photon spectral index for Mrk 766 (ax =- 
1.5) in their sample of NLSl galaxies, yet steeper tha n what most 
broad line Seyfert 1 or even some other NLS 1 show. lOrupe et alj 
i200ll) give a ax =-1.77 for this source, which is exactly the mean 
value for the whole sample of soft X-ray selected AGN reported in 
[Grupe 1 2004), not as flat as initially found. Nandra et al. ( 1997) re- 
ported the detection of an Fe Ka line in the X-ray spectrum. More- 
over, recent XMM-Newton observations have caused a debate over 
whether or not the soft X-ray spectru m arises from a du sty warm 
absorber or relativistic emission lines ^M ason et all2003h . 

In the ultraviolet, the spectrum of Mrk 766 looks also un- 
usual among Seyfert 1 galaxies due to the weakness of both th e 
ultraviolet continuum and big blue bump <Walter & Finkl ll993h . 
The deficit of ultra violet photons has been attributed to redden- 
ing (Page et al. 1999) and/or to absorption I Walter & Fink''l993'). 
Actually, Walter & Fink ( 1993) claim that there is no basic differ- 
ence between the unabsorbed optical to X-ray continua emitted by 
Mrk 766 and by other Seyfert I galaxies. The dusty scenario is 
supported by the optical polarization of Mrk 766, of ~2 per cent 
increasing to the blue, due to scattering by dust grains. The grains 
are probably located within the narrow line region giv en that the 
broad lines show more polarization than the narrow lines iGoodricll 
Il989l) . 

In the near-infrared, much of t he information available for 
Mrk 766 comes from imaging studies. lAlonso-Herrero et ^il998h 
report that the surface brightness of this source at the 0.5-1 kpc 
radius is I-2magarcsec~^ brighter than that of normal galaxies. 
Moreover, it shows one of the reddest H-K colours for the 0.5 and 
I kpc diameter apertures of their sample. This result is confirme d 
by the J- and A"-band imaging work of iMarguez et alj Jl999h . 
where it is shown that Mrk 766 displays a double nuclear struc- 
ture with a (J-K') color index redder than the rest of the galaxy. 
That structure seems to correspond to a feature delineated by the 



dust pattern that surrounds the innermost 2" in the HST optical 
image (Malkan et al. 1998). Spectroscopically, Mrk 766 have been 
poorly studied in this band. To our knowledge, the only observa- 
tions of Mrk 766 published in the NIR region is the A^-band spec- 
troscopy of Grupe & Thomas (2002) and the L-band spectroscopy 
of iRodrfguez-Ardila & Viegas. Ii2003>) . In the former, the spectrum 
clearly reveals Br7 and probably the H2 and [Si VI] lines. In the 
later, the 3.3/^m emission, attributed to polycyclic aromatic hydro- 
carbons (PAH), is detected. The observation of PAH emission is 
usually taken as unambiguous evidence of starburst activity. Inter- 
estingly, a strong starburst might be behind the high [O IIl]/H/3 ra- 
tio observed in Mrk 766, so extreme (Grupe 2004) as to possibly 
exclude M rk 766 from the NLS 1 class, following the definition of 
NLSls bv iOsterbrock & Pogge (.1985 ). 

For all the above, Mrk 766 arises as a very interesting target 
to continue our analysis of the continuum and NLR properties of 
NLSl galaxies. To this purpose, NIR observations supplemented 
with archival UV - optical and mid-far infrared data will be used in 
this work. This set of data is described in Sec.|2| The description 
and characterization of the observed NIR continuum is described 
in Sec.|3| A discussion about the most important emission line fea- 
tures detected in the spectrum, both from the BLR and NLR, is 
presented in Sec.|4] In a first approach to apprehend the nature of 
Mrk 766, we compare its SED with that of NLS 1 galaxies (Sec.|5j, 
and, particularly, with the SED of Ark 564, which has been al- 
ready modeled (Paper I ). This step will help at constraining the 
models for the evaluation of the emission-line ratios presented in 
Sec.|6| After modeling the emission-line spectrum, the associated 
(and consistent) SED of the continuum is discussed in Sec.0 The 
concluding remarks appear in Sec.|8| 



2 OBSERVATIONS 

NIR spectra of Mrk 766 in the interval 0.8—2.4 ^m were ob- 
tained at the NASA 3 m Infrared Telescope Facility (IRTF) on the 
nights of April 21 a nd 25, 2002 (UT) using the SpeX spectrometer 
teavner et al]|2003h . In addition, spectra covering the interval 3- 
4 jj,m were also obtained with the same setup on the night of April 
22 (UT). In both cases, the detector consisted of a 1024x 1024 AL- 
ADDIN 3 InSb array with a spatial scale of 0.15"/pixel. Simulta- 
neous wavelength coverage was obtained by means of prism cross- 
dispersers. A 0".8x 15" sht, oriented at the PA of 112°, was used 
during the observations. The resulting spectral resolution was 360 
kms~^ for both setups. The observations were done nodding in 
an ABBA source pattern with individual integration times of 120 s 
and 30 s for the short (0.8-2.4 ^m) and long (3^ ^im) wavelength 
range, respectively. The total on-source integration time was 1680 s 
and 1080 s for the April 21 and 25, respectively, and 900 s for the 
22 April observation. The AO V star HD79108 (V=6.13) was ob- 
served after the target's to remove the telluric bands. It was also 
used to flux calibrate the spectra. The spectra were reduced fol- 
lowing standard procedures for NIR spectroscopy. They include 
subtracting consecutive 2-D AB pairs to remove background, flat 
fielding the resulting imaging, extracting the signal along the dis- 
persion axis, co-adding and calibrating. The signal within the cen- 
tral 0".8x 1" was summed up to obtain the nuclear spectrum. The 
1" extraction window along the slit, with centre at the peak contin- 
uum luminosity in K, is in accord to the seeing of 0."8 measured 
during the observations. This window size corresponds to the in- 
ner 250 pc (Ho=75 km s'^ Mpc"^) of the AGN. No effort to ex- 
tract information from the extended emission, seen up to ~400 pc 
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Figure 1. Observed NIR spectium of Mrk 766 in laboratory wavelengths. In order to easy the identification of the spectral features, the most important 
permitted and forbidden lines are marked. 



NW and SE of the centre, was done. The reduction and extrac- 
tion procedures were done using the in-house software Spextool 
iCushine et al.lE()04l) provided by the IRTF observatory. For tel- 
luric band correction and flux calibration, a script also provided by 
IRTF, designed to this purpose (Vaccaet al. 2003) was employed. 
Afterward, the two spectra in the short wavelength range were av- 
eraged to form a single 0.8—2.4 spectrum. 

FigureQshow the final rest-frame NIR spectrum of Mrk 766. 
Note that within the intervals 1.74 ^m-1.86 /im and 2.4 /im- 
2.9 /im, the atmospheric transmission drops to zero. For this rea- 
son, these two regions are missing in FigureQ As can be seen, the 
spectrum is dominated by permitted emission lines of H I, He I, O I 
and Fe 11. Also prominent are the forbidden lines of [S III] AA9068, 
9531, [Si VI] 1.963 fim, [Mg VIII] 3.03 fim and [Si IX] 3.936 ^m. 
Note also the bum p centred at 3.3 nm, attributed to emission of 
PAH molecules I R odriguez- Ardila & Viegasll2003h . It represents 
an unambiguous signature of starburst activity within the central 
250 pc of this NLS 1 galaxy. 



Fluxes of the identified NIR lines measured from our spectrum 
are listed in TableQ They were calculated by fitting a Gaussian to 
the observed profile a nd then integrating th e flux under the curve. 
The LINER software jPogge & Owenll993h was used for this pur- 
pose. 

Archival HST data taken with the Space Telescope Imaging 
Spectrograph (STIS) on 1 1 April 2000 (UT) were also used in or- 
der to measure the flux of the most important UV lines of Mrk 766. 
The spectra were obtained through the 52" x 0.2" slit and the low- 
resolution (R~750) G140L and G230L gratings. The spectra cover 
the rest-frame wavelength range 1100-1690A and 1550-3200A, 
respectively. The final rectified wavelength calibrated x2d spectrum 
produced by the //5r pipeline was used to extract the 1-D spectrum. 
The extraction window in the spatial direction was the same used 
in the NIR data. The IRAF package was used to this purpose. 

In the optical region, fluxes for the most i mportant lines were 
taken from Gonzalez Delgado & Pered il996h (hereafter, GP96). 
They report relatively strong Fe II emission and coronal lines, as 
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well as extended optical emission through a region the total size of 
8" (about 2 kpc) along the stellar bar. Because of the high quality 
of the data and resolution, that matches ours', the nuclear fluxes of 
GP96 will be used through this work. 

Finally, the fluxes of the emission lines from UV to NIR 
were derredened for Galactic and internal extinction. The v alue 
of E( B-V) corresponding to our Galaxy is 0.02 ( Schleeel e t alj 
Il998t) . The internal extinction E(B-V)=0.5±0.1, applied to the 
data, was determined from the average of the extinction values de- 
rived from the flux ratio of the narrow components of Pa/3/Pa7, 
E(B-V)=0.69, and that determined by GP96 [E(B-V)=0.38]. The 
value found from the Pashen lines agre es with that derived from 
the Balmer decrement iGrupe e t al. 1998*). indicating that Mrk 766 
is indeed highly reddened in the optical/UV. Note that higher val- 
ues of extinction were also found, for example, from the flux ra- 
tios [Fen] 1.257^m/1.644^m [E(B-V)=1.2] and Br7/Pa/3 [E(B- 
V)=1.15]. However, we rejected the latter two values because it 
is not certain if the bulk of the NLR gas is affected by such large 
values of extinction, mainly if we take into account that the optical 
spectrum of Mrk 766 displays prominent [Ne V] 3346, 3426 Alines 
(see Figure 3 of GP96), which are located in the blue end of the 
spectrum. Table|2|list the dereddened o ptical and NIR na rrow line 
fluxes. Note that the extinction law of lCardelli et aljri989ll) was ap- 
plied in that correction. 

The choice of the lCardelli et alJ ^98^ extinction law is due 
to the fact that we are assuming that the dust affecting the NLR is 
located in a non-nuclear environment. It is yet possible that the ob- 
served continuum in the UV-optical region is affected by the same 
amount of extinction, but if there is an additional nuclear dust com- 
ponent, it is likely that the Cardelli et al. ( 198^ law is no longer 
valid. The main support to this is the lack of the A2175 dip in 
the UV conti nuum shape. This is in accord to recent claims by 
iGaskell et alJ (2004), who provided substantial evidence of a rel- 
ative lack of small grains in the nuclear dust, which would explain 
the absence of the A2175 dip in most AGN spectra. They proposed 
a reddening curve fairly universal, being quite flat in the UV but se- 
lective in the optical. Due to the lack of enough continuum points in 
the optical region that adequately connects the NIR with the HST 
data, we feel that any correction to the nuclear UV continuum is 
highly uncertain. Since our interest is concentrated in the NLR gas 
and the reprocessed continuum by this gas, we feel that this correc- 
tion does not affect our main results and conclusions. 



3 THE NIR CONTINUUM 
3.1 The observed continuum 

The continuum emission displayed by Mrk 766 in the interval 
0.8—4 /.im is complex (see Figure|2j. At least three main features 
can be easily identified in its spectrum: (i) a steep continuum blue- 
ward of 1 fj,m, very similar to the the so-called small blue bump in 
the optical region observed in AGN around 4000 A; (ii) a turnover 
of the continuum at ~ 1.1 ^m. The continuum emission raises long- 
ward up to ~2.2 /^m, where it reaches a maximum; (iii) a steep de- 
cline of the continuum flux with wavelength, longward of 2.2 ^m. 

Few AGN at a similar wavelength coverage and spectral res- 
olution in the NIR region are available in the literature, so a 
meaningful comparison of the continuum observed in Mrk 766 
with that of other AGN is not pos sible. However, data pub- 
lished for the NLS Is galaxies IZw 1 jRudv et aljr2000l) . Mrk 478 
jRudv et alJ boOlh and Ark 564, Ton S 180, IH 1934-063 and 
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Figure 2. (a) NIR SED observed in Mrk 766. The dotted line is best-fitted 
power-law function. Note the strong excess of emission over the power- 
law, in the interval 12000-32000 A. Note also that the power-law overpre- 
dicts the continuum emission lonwai-d of 38000 A. (b) NIR SED modeled 
with a composite fu nction: a power-law (dotted lines), of spectral index 
o=-0.44, as found bv lAlonso-Herrero et alj|2003l) . plus a blackbody distri- 
bution (dashed line) of T\,i,=\200 K. The thick line represent the sum of 
both components. 



Mrk 1044 iRodrfguez-Ardila et ai]l2002ei) show that the turnover 
of the continuum near I.I/xm is a common feature and may repre- 
sent the shift from the non-thermal continuum emitted by the active 
nucleus to the stellar and thermal dust emission from the circumnu- 
clear region that dominates at longer wavelengths. The rapid rise of 
the continuum in the I and H bands with a peak at the K-band, and 
then its slow decline, has not been reported before, to our knowl- 
edge, f or an AGN but I Zw 1 . Note, however, that the I Zw 1 spec- 
trum of lRudv et ail j200(il ends at 2.40 ^m. 

In order to characterize the NIR continuum emission in 
Mrk 766, a power-law of the form F^ oc A" was fitted to the data. 
The task nfitld of the STSDAS package of IRAF was used to this 
purpose. Special care was taken to not include emission lines in 
the spectral windows used in the fit. The result is shown in the 
upper panel of Figure |2| The derived spectral index a was -0.16, 
significantly fla tter than the average valu e of -0.64 found for the 
AGN sample of lEdelson & Malkajh986l) and also fl atter than the 
value of -0.44 found bv lAlons^Se^iroet alj i2003h in the inter- 
val 1 — I6^m for Mrk 766. It can be seen from Fig.|2|that the fitted 
function overpredicts the continuum emission redward of 4/im. In 
addition, it shows that a single power-law cannot describe the con- 
tinuum emission in the 1— 4^m region. There is an excess of emis- 
sion over the power-law that in no way follows a simple function. 

A closer look to the overall shape of the excess of emission 
over the power-law, in the wavelength range 1-4 fim, suggests that 
it approaches to that of a blackbody distribution. In order to test 
this hypothesis, we fitted a composite function (power-law plus 
blackbody) to the observed NIR SED. In the fit, the spectral index 
of the power-law was constrained to the value q:=-0.44, found by 
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Table 1. Observed NIR emission Une fluxes measured in Mrk766 





Flux 




Flux 


Line 


xlQ— "'^^ erg cm" ^ 


Line 


xlO""*^^ erg cm 


Bra 4.052/im 


31.3±4.70 


Hell(b*) 1.162/im 


2.94±0.81 


[SilX] 3.935/im 


17.06±3.60 


[PlI] 1.146/im 


1.44±0.29 


[MgVIII] 3.027^m 


26.6±8.06 


Ol 1.128^m 


48.86±2.40 


[CaVIII] 2.321/im 


2.47±0.6() 


Pa7(n*) 1.937^tm 


13.05±0.77 


H2 2.223^tm 


0.81±0.15 


Pa7(b*) 1.937(Um 


87.30±3.07 


Br7 (n*) 2.157^tm 


7.4()±().35 


Fell 1.086/im 


50.6±2.40 


Br7(b*)2.157/im 


20.0±2.26 


Hel 1.083/:/m 


216.4±1.53 


H2 2.121/im 


2.36±0.24 


Fell 1.050/im 


14.86±1.96 


Hel(n*) 2.058Mm 


2.03±0.20 


[So] 1.037/im 


1.67±0.46 


Hel(b*) 2.058Atm 


6.09±1.10 


[Sll] 1.033Mm 


2.11 ±0.46 


H2 2.033^tm 


1.73±0.22 


[Sll] 1.032/im 


3.38±0.46 


[SiVI] 1.963/im 


15.6±0.94 


[Sll] 1.028/im 


2.91±0.46 


H2 1.957^tm 


4.39±0.39 


Fell l.()17^m 


6.23±2.10 


[Fell] 1.634/im 


8.20±0.40 


Hell(n*) 1.012/im 


18.57±2.02 


[Six] 1.430Atm 


6.47±0.54 


Hell(b*) 1.012/im 


5.68±0.70 


Pa/3(n*) 1.281/im 


27.15±0.50 


Pa5(n*) 1.004/im 


4.82±0.69 


Pa/3(b*) 1.281/im 


117.8±1.80 


Pa5(b*) 1.004/im 


68.7±2.7 


[Fell] 1.257/im 


7.69±0.70 


Fell 0.999/im 


21.98±2.13 


[Six] 1.252/im 


9.42±0.96 


[Svill] 0.99lAtm 


lO.Oitl.OO 


[Pn] 1.188Atm 


2.19±0.29 


[Cal] 0.985^jm 


5.90±1.34 


Hell(n*) 1.162/im 


1.40±0.24 


[SlII] 0.953Atm 


110.0±0.90 



n and b stand for the narrow and broad components, respectively, of the permitted emission hnes. 



Table 2. Line fluxes in Mrk766 corrected for E(B-V)=0.5 and normalized to F6007+ x 100 



Line 


Hux 

10-" ergcm-2 s"! 




Line 


Hux 

10-" ergcm-2 s'^ 


Fa^5007 


[Silx] 3.935Mm 


1.87±0.39 


1.10±0.25 


[SlI] 6717 A 


4.14±0.33 


2.45±0.26 


[MgVIII] 3.027/im 


2.94±0.88 


1.74±0.54 


Hel 6678 A 


2.00±0.22 


1.18±0.16 


[CaVIII] 2.321/im 


0.30±0.07 


0.17±0.04 


[Nil] 6584 A 


18.01±1.08 


10.66±0.99 


Hel 2.058Mm 


0.96±0.13 


0.57±0.09 


Ha 6563 A 


281.1±16.9 


166.3±15.5 


[SiVI] 1.963/im 


1.87±0.11 


l.lliO.lO 


[Nil] 6548 A 


4.00±0.12 


2.37±0.18 


[Fell] 1.634/im 


1.05it0.05 


0.62±0.05 


Fex] 6375 A 


0.45: 


0.27±0.02 


[Six] 1.430;um 


0.91±0.07 


0.54±0.06 


[Oi] 6363 A 


1.28±0.12 


0.76±0.09 


Pa/B 1.281/im 


21.34±0.21 


12.63±0.91 


[Oi] 6300 A 


2.80±0.22 


1.66±0.18 


[Fell] 1.257 urn 


1.12±0.10 


0.66±0.08 


[Fevii] 6086 A 


2.40±0.12 


1.42±0.12 


[Six] 1.252/im 


1.33±0.13 


0.79±0.10 


Hei 5876 A 


13.74±1.10 


8.13±0.87 


[Pii] 1.188/tm 


0.35±0.05 


0.21±0.03 


[Nil] 5755 A 


0.37±0.06 


0.22±0.04 


Hell 1.162/^m 


0.68±0.13 


0.40±0.08 


[Fevil] 5721 A 


1.33±0.12 


0.79±0.09 


[PlI] 1.146^jm 


0.23±0.05 


0.13±0.03 


[Ol] 5577 A 


0.33±0.07 


0.20±0.04 


Pa7 1.937^m 


12.27±0.37 


7.26±0.56 


[Olll] 5007 A 


169.3±11.9 


100.0±7.10 


Hel 1.083//m 


36.02±0.36 


21.31±1.53 


H/3 4861 A 


82.93±2.95 


49.±4. 


[Sll] 1.037/im 


0.29±0.08 


0.17±0.05 


[Arvi] 4740 A 


0.67±0.09 


0.39±0.06 


[Sll] 1.033//m 


0.36±0.08 


0.21±0.05 


[Arvi]4711 A 


0.63±0.09 


0.37±0.06 


[Sll] 1.032/nm 


0.57±0.08 


0.34±0.05 


Hen 4686 A 


7.34±0.51 


4.34±0.43 


[Sll] 1.028/Lim 


0.53±0.09 


0.32±0.06 


[Oiii]4363 A 


3.47d=0.59 


2.05±0.38 


Hen 1.012/im 


4.27±0.38 


2.53±0.29 


H7 4341 A 


34.82±3.13 


20.6±2.36 


Pa<5 1.004^m 


12.94±0.52 


7.66±0.62 


[SlI] 4070 A 


2.54±0.46 


1.50±0.29 


[SviII]0.991^m 


1.87±0.19 


1.11±0.14 


[NeIII]+Hl 3968 A 


6.27±0.63 


3.71 ±0.46 


[Cal] 0.985^m 


1.04±0.24 


0.62±0.15 


H8+Hel 3888 A 


3.60±0.61 


2.13±0.39 


[Slll] 0.953/im 


20.54±0.21 


12.16±0.87 


[Nelll] 3868 A 


20.41 ±1.63 


12.08±1.29 


[Nill] 7378 A 


0.40±0.04 


0.24±0.03 


H9+Heii 3835 A 


1.09±0.51 


0.65±0.31 


[On] 7320+7330 A 


1.87±0.11 


l.llitO.lO 


[Fevii] 3758 A 


4.00±0.92 


2.37± 0.57 


[Arm] 7135 A 


3.20±0.22 


1.89±0.19 


[On] 3727 A 


22.28±1.56 


13.18±1.31 


Hel 7065 A 


3.07±0.21 


1.82±0.18 


[Nev] 3426 A 


28.01±1.68 


16.58±1.54 


[Arvi] 7006 A 


1.60±0.11 


0.95±0.09 


[Nev] 3345 A 


10.0±0.80 


5.92±0.63 


[Sll] 6731 A 


5.34±0.32 


3.16±0.29 
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lAlonso-Herrero et alJ jlOO^) while the temperature and amplitude 
of the blackbody function were left as free parameters. The result, 
displayed in the lower panel of Figure |2| show that the composite 
function, with a blackbody of temperature rbb=1200 K, adequately 
describes the continuum distribution. Integrating the flux over these 
two components in the observed interval gives Fbb=l-55x 10~^^ 
ergcm"^ s^^ and FpL=5.60x 10~^^ ergcm"^ s~^ for the black- 
body and power-law functions, respectively. It means that in the 
NIR region, the former component accounts for 28% of the total 
continuum flux. 

What is the origin of each of the two continuum components, 
both emitted in the inner 250 pc of the AGN? Given the composite 
nature of Mrk 766, harboring both a dusty active nucleus and a cir- 
cumnuclear starburst, there are several possibilities. The power-law 
continuum can be associated to continuum emission from the AGN 
plus some contribution from circunnuclear stellar population and 
dust. Note that this power-law cannot simply be regarded as an ex- 
trapolation to the power-law used to describe the continuum in the 
UV and optical region. This is because in the optical, the continuum 
displayed by Mrk 766 is essentially flat, even in the UV region (see 
Figure 3 of GD96). This suggest that the continuum emission, at 
shorter wavelengths, is strongly reddened. Under the assumption 
that the power-law found in the NIR carries a significant contribu- 
tion from the AGN, it can be used to find out an upper limit to the 
extinction affecting the UV and optical region. We determined that 
an at least a E(B-V)=0.3 to the continuum is necessary in order to a 
power-law of spectral index a=-0.44 can adequately represents the 
AGN continuum in the optical region. 

Regarding the nature of the excess of emission well de- 
scribed by a black body fu nction, it could be relate d to a compact 
dust/molecular thick torus iGranato & Danese* '1994') or result from 
emission from hot dust (T>1000 K). This reasoning follows sim- 
ilar distributions found in other AGN by iMarco & AUoirJ IT99I 
and lMarco & AUoiiJ <2000l) . 



3.2 The hot dust hypothesis 

The excess of emission over the power-law distribution in the 
NIR region and its adequate representation obtained after adding 
a blackbody distribution suggest that it can be due to hot dust 
grains. The temperature derived from the fit (T'bb=1200 K) is near 
to the peak value for the hottest possible grains, i.e., graphite 
grains at their evaporation temperature, T~1500K, and higher 
than the sublimation tem perature of silicate grains (T~IOOO K; 
lOranato & Danes3ll994^ . Considering that our spatial resolution 
is limited to ~250 pc, it is very likely that dust at higher temper- 
atures exists closer to the central source, ruling out the possibil- 
ity of silicates as the main component of the nuclear dust grains. 
This result does not exclude that silicate grains can exist within the 
NLR of Mrk 766. They still can be present within the observed 
region but a much lower temperature. Recalling that HST observa- 
tions of Mrk 766 reveals filaments, wisps and irregular dust lanes 
around the unresolved nucleus ( iMalkan et all 19981) and that the op- 
tical polarization of Mrk 766 is due to scattering from dust grains 
jGoodriclJll989h . the dusty NLR scenario for Mrk 766 is strongly 
supported by the NIR data. 

Using the temperature of the blackbody distribution derived 
above as the average temperature of the graphite grains and a K 
band flux of 1.03 x 10"^^ erg s"^ cm"^ Hz"^ at 2.2 /^m found for 
the blackbody component after subtracting the power-law, we can 
roughly estimate the dust mass associated with the excess of emis- 



sion. Following |Barvaiiii3 ^983), the infrared spectral luminosity, 
in ergs s~^ Hz~^, of an individual graphite grain is: 

Ll'i, = ^TT nQ.B^T^r) (1) 

where a is the grain radius, Q^. — qirv'' is the absorption 
efficiency of the grains and Bv{Tgi) is the Planck function for 
a grain of temperature Tgr. Adopting, as in|Barvainis (1987), a 
value of a=0.05/im for graphite grains and Q^=0.058 and setting 
rgr=1200 K, we find Lf .=9.29x 10"^* ergs s"^ Hz~\ 

The total number of emitting grains (hot dust) can be approx- 
imated as: 

JVhD ~ ,gr ■ (2) 

Finally, for graphite grains, with density pg=2.26 g cm~^, 
47r •< 

Mhd ^ — a'^ATHDPg (3) 

Taking Mrk 766 at a distance of 51.6 Mpc (z=0.0129, 
Ho=75 kms"VMpc), we obtained AfHi3=3.54x lO'*''' and 
Mhi3=2.1x10-=' M0. 

The mass of hot dust derived fo r Mrk 766 i s similar to that 
found in NGC 37 83 (2.5 xlO'^ MQ.lGlassll992h and NGC 1068 
(1.1x10"^ M0,Im arco & Alloin"200Q'); larger than that derived 
for NGC 1566 (7x10"'' M0, Baribaud et al. 1992) and NGC 4593 
(5x10""' M0, Sant os-Lleo et al . 1995) and significantly smaller 
than that of NGC 7469 (0.05 M0, Marco & AHoin 1998) and 
Fairall 9 (0.09 M0. lciavel et alJl989ft yNote that for all the above 
objects but Mrk 766, the mass of hot dust was derived from photo- 
metric observations. The spectral resolution of our data offers the 
possibility of examining a spectral feature not previously observed 
in such a detail. 

Regarding to the location of the hot dust, our spectral reso- 
lution imply that it must reside within the inner 250 pc from the 
centre, since this is the region size covered by the nuclear NIR 
spectrum employed in this work (see §|2j. The high temperature 
found for the dust, however, allows us to impose a tighter con- 
straint to the location of the emitting region: it would be located 
in the central 100 pc. This value is deduced by inspecting Fig. 3 
of Marco & Alloin 1 1998), where a plot of dust temperature as 
function of distance from the central source was constructed for 
NGC 7469. Note that in this reasoning, we have assumed that the 
dust in Mrk 766 has similar properties to that of NGC 7469. Fur- 
ther support to this value can be obtained following the results for 
NGC 1068 I Marc o & All oin 2000). These authors concluded that 
hot dust (rgr=I500 K) should be extremely confined and located 
at a radius less than 4 pc. If that is also the case for Mrk 766, it 
is most likely that the excess of emission in the NIR is associated 
to the torus than to circumnuclear dust. Only high resolution spec- 
troscopy will be able to disentangle this possibility. Whatever the 
case, our results provides the first direct spectroscopic evidence of 
hot dust in AGN and show the potential that NIR spectroscopy has 
at unveiling that component. 



4 THE LINE EMITTING REGIONS 

A careful inspection to the spectrum of Mrk 766 shown in Fig. Q 
evidences prominent permitted emission lines of Fe II, Hi, He I 
and Ol and bright forbidden lines of [Slll]. These features carry 
most of the energy radiated in the form of emission lines in 
the 0.8 — 1.3/im interval. In addition, several coronal lines with 
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Figure 3. Top: Mrk766 spectrum with the most important emission lines identified in the region between 0.8—1.1 /im. The continuum emission has been 
subtracted and the resultant spectrum displaced by an arbitrary amount for displaying purposes. The thick line is the best empirical Fe II template derived 
from the models of Siaut & Pradhan 1 2003 ). Note the good agreement between the template and most of Fe II features, mainly the 9200A blend and the 1 fim 
lines. Bottom: Fe II Empirical template with the most important lines identified. Nothe that in order to ease the identification of individual lines, we show the 
template before being convolved with a Gaussian that matches the width of individual permitted lines of Mrk 766. 



ionization potential (IP) of up to 350 eV, some of them previ- 
ously observed in just a handful of AGN, are detected. Such is 
the case of [Svill] 0.991 /xm (1P=280.9 eV), [Six] 1.252 fim 
(IP=328.3 eV), [Sivi] 1.963 fim (IP=166.7 eV), [Siix] 3.936 /xm 
(IP=303.17 eV), [Six] 1.143 fim (IP=351.1 eV) and [MgVIIl] 
3.03 (IP=224.95). Low ionization hnes, with 1P<13.6 eV, are 
also present (i.e., [C I], [Fe II] and [S II]). 

The simultaneous observation of very high and neutral/low 
ionization lines imply that a wide variety of physical conditions 
must coexist within the central few hundred parsecs of this object. 
In the following sections, with the help of ionization models, we 
will extract information about the structure and physical state of 
the NLR, taking as reference, the observed continuum and emis- 
sion line spectrum. A detailed analysis of the physical conditions 
and modeling of the BLR features will be discussed in a separate 
paper (Rodriguez-Ardila et al. 2004, in preparation). Here, we will 
briefly mention the most important features identified in the BLR 
spectrum of Mrk 766. 

Few AGN with positive d etection of N IR Fe II lin es have been 
reported in the literature so far iRudv etal. 2000, 20o'll: lNagar et alJ 
l2002t iRodriguez-Ardila et aljl2002'al) . Here, we report the identi- 



fication of the Fe II 9200 A blend and the Fe II 1 /im lines^ in 
Mrk 766, as is shown in Fig. |3| In the upper part of the Figure, 
we have marked the strongest transitions identified. Overplotted, 
is the empirical template we derived. This work was done with 
the he lp of the NIR Fe II line list and models of Sigut & Pradha^ 
(2003). A zoom of the empirical template, plotted at the bottom, 
was indeed const ructed using the positio n and relative intensity of 
the lines listed bv lSigut & Pradhajj2003h . We then adjusted by eye 
most of the line intensities in disagreement with the model in order 
to get a better approximation to the observed spectrum. 

The detection of NIR Fe II lines provides important informa- 
tion on the mechanism driving the formation of this ion. Primary 
cascade lines descending from the upper 5p levels to the lower 
e*D and e^D terms, thought to be pumped by hya fluorescence 
(Johansson & Jordan 1984; Sigut & Pradhan 1998, 2003), are lo- 
cated in the region between 8500 A and 9500 A. Their detection, as 
in Mrk 766, supports this mechanism as a competing one in the for- 
mation of the Fe II. Moreover, the l^m lines can also be pu mped by 
Lya fluorescence (see Figure 16 o f ISigut & Pradhanl2003h . The en- 



^ It is known under this name the set of Fell lines located between l/im 
and l.lfim, among which 0.9997fj,m, 1.0I71/im, 1.050/.tm, 1.0863/im and 
1.1 126/im are the strongest ones 
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Table 3. FWHM, in km s'l, of the detected lines in Mrk766 



Line 


FWHM 


Line 


FWHM 


LV^d. VlllJ i. fJilll 




ne 11 i^ii^ i.io^iii 








ne 11 \u) i.io^iii 


1220 


Dl ^ \LL J 




rp ttI 1 1 A/irrt 

[r iij i-it^iii 






1600 


T 1 1 ^R/zm 
\J 1 l.lZo/Jilll 


1270 


H2 2 212 


^360 


Pa'V ^^n^ 1 00'^(;m 


^360 


He I (n) 2.058^m 


<360 


Pa7(b) 1.093^™ 


1720 


He I (b) 2.058^m 


1890 


He I 1.083/tm 


780 


H2 2.033//m 


<360 


[S II] 1.037/im 


^360 


[Si VI] 1.96^m 


390 


[S II] 1.033^im 


^360 


H2 1.957^m 


^360 


[S II] 1.032^tm 


^360 


[Fe II] 1.63//m 


<360 


[S II] 1.028/im 


^360 


[Si X] lA3fim 


430 


Hell(n) 1.012/im 


^360 


Pa/3(n) 1.281/im 


<360 


He 11(b) 1.012/im 


1250 


Pa/3(b) 1.281/im 


1740 


Pa<5 (n) 1.004/im 


^360 


[Fe II] 1.257^m 


370 


Pa(5(b) 1.004/im 


1700 


[S IX] 1.252/im 


600 


[S VIII] 0.991^tm 


^360 


[PlI] 1.18/zm 


^360 


[Ca I] 0.985/tm 


^360 






[S III] 0.953^tm 


^360 



ergy levels that are populated after their emission originate nearly 
50% of the optical Fe 11. This de-excitation channel has been over- 
looked for years and should be included in modeling the Fe II spec- 
trum. Given the importance that the Fe II lines have in AGN studies, 
a paper describing the detection of NIR transitions and the physi- 
cal information that can be derived for a sample of AGN, including 
Mrk 766, is currently in preparation. 

In regard to the NLR, the width of the lines provide impor- 
tant clues about the velocity field of the emitting gas. This infor- 
mation will be useful for the models that are discussed in Sees. |5| 
and|6| Table|3|lists the FWHM measured for the lines detected in 
the spectrum. 

The values of FWHM found in Mrk 766 indicate that all 
NLR lines but [Six], [Si VI] and [Six] are spectroscopically un- 
resolved or in the limit of the spectral resolution. This agrees with 
the FWHM of the optical lines measured by GDP96. They found 
that H/3, [OIII], Ha, [Nil], [S II] and [S III] had widths of less 
than250kms~^,wellbel ow our spectral resolution of 360 km . 
IVeron-Cettv et alj <200lh also reported a FWHM of 330 km s"^ 
for [Olll] A5007. This supports the existence of two distinct re- 
gions within the NLR of Mrk 766: one in which low and medium 
ionization lines are formed (the classical NLR), and another where 
the bulk of coronal lines (the so-called coronal line region, CLR) 
are emitted. Note, however, that no asymmetries are observed in the 
line profiles of the NIR coronal lines, consistent also with the obser- 
vations of GDP96. This contr ast to what is usually found in high- 
ionization lines in AGN (.Penston et alJll984 lErkens et alJll997t 
iRodrfguez-Ardila et alj|2()04bl) . The lack of asymmetries suggest 
that the coronal lines are weakly related to outflowing material 
or strong shocks, as is usually claimed lErkens et al. 1997). They 
should mostly originate from interaction between matter and radia- 
tion from the central source. This issue will be discussed in Sec.|6| 

In addition to the permitted and forbidden lines, the NIR spec- 
trum of Mrk 766 also displays conspicous molecular H2 lines, par- 
ticularly in the K-ban d. Analysis of this emission is presented in 
iRodrf guez- Ardila et al . (2004a), where it is concluded from the 
analysis of molecular line ratios and the vibrational and rotational 
temperature that thermal mechanisms drive the excitation of the 
H2 gas. Moreover, X-rays from the central source can account for 




log(i/) 



Figure 4. Comparison of the observed spectral energy distribution of Mrk 
766 (open circles) and the model of Ark 564. The curves correspon d to the 
fit to th e Ark 564 continuum obtained by detailed modeling 1 Co ntini et alJ 
l2003tJ) . The dash-dotted line represents the old stellar population back- 
ground. The vertical dotted lines refer to the IR wavelengths at 10, 16, 25, 
60 and 100 ^im. 

nearly 70% of the observed emission. The remaining 30% was at- 
tributted to the starburts component, mostly via shock heating by 
supernova remnants. The molecular spectrum will no be further 
discussed throught this text. 



5 MODELING OF THE BROAD-BAND SPECTRAL 
ENERGY DISTRIBUTIONS 

In order to model the SED of Mrk 766, some indication on which 
type of clouds may be present in the nuclear region is required. A 
comparative analysis of the observed continuum of Mrk 766 to that 
of other AGN provides a first approach to this problem. 

The observed continuum spectra of Mrk 766 (open circles) is 
plotted in Fig.|4| As a guide, the energy distribution obtained in the 
modeling of Ark 546 (Paper I) is also plotted. The vertical dotted 
lines refer to wavelengths in the IR at 10, 16, 25, 60, and 100 fj,m. 

The d ata come from the NED^. Observations in the optical- 
UV range (G ioia et alJll99(tlGavazzi & Bosellilll996llGioia et"!?] 
1 199(1 iDe V aucouleur set alJll99lh have been corrected for inter- 



nal reddening (E(B-V)=0.5). The other data refer to observations at 


lower frequencies (Balzano & Weedman 1981; 


Dressel & Condon 


1978; Moshir et al. 1990; Becker et al. 1995; Maiohno et aljl993. 


iRudv 


et aL 


■ 1982bs; 


Soinoglio et aL 


■ 1995; .Alonso-Herrero et alJ 


I2OO3I 


. Actually, Alonso-Herrero et al. claim that the nuclear fluxes 



at wavelengths greater than approximately 2 ^m will not be signif- 
icantly affected by moderate amounts of extinction. Moreover, the 
view of the torus through a fragmented screen of dust will dominate 



^ Nasa Extragalactic Database, |http://nedwwww.ipac.caltech.edu| 
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the flux in the near-IR and simulate a low level extinction. No data 
appear in the UV where the corre ction could be signifi cant. 

The X-ray data, taken from 'M ason et alJ l200?i) using XMM 
observations, ha ve been correcte d for absorption using the coeffi- 
cients given by jZombecklll99d) adopting a column density Nh = 
10^^ cm-^ 

In Fig. |4] the Mrk 766 SED is compared to the continuum 
spectrum obtained from the detailed modeling of Ark 564 (Paper 
I). Recall that the continuum which results from modeling rep- 
resents reprocessed radiation (free-free and free-bound radiation) 
from clouds reached by the flux from the central source and heated 
by the shocks. The main reason for this comparison is the fact that, 
overa ll, both objects have similar spectral properties ( Page et alj 
1 19990 . For example, they do not show evidence for the Big Blue 
Bump between the ultraviolet and the soft X-ray region. Moreover, 
they show a low ratio between the ultraviolet and the soft X-ray 
fluxes, implying a fiat optical - X-ray energy distribution, as eas- 
ily seen in Figure|4] In fact, for Mr k 766 and Ark 56 4, aox equals 
to -1.17 (tiiis paper) and -1.08 (.Romano et aljE004l) . respectively, 
indeed flatter than the mean va lue of -1.63 for low-redshift radio- 
quiet objects jYuanetalll 99i). It suggest that, in a first approach, 
the ionizing sources of Ark 564 and Mrk 766 are very similar 

Note, however, that the flux measurement in the X-rays re- 
gion of Mrk 766 at the highest frequencies are lower than that for 
Ark 564. This may indicate that shock velocities as high as found in 
the NLR of the later (V^ 1000 km s~^) are likely not to be present 
in the former. The data do not exclude a contribution to the X-rays 
coming from clouds with Vs of about 500 kms^^. These type of 
cloud s are found to be present in Seyfert 2 galaxies iContini et alj 
l2004h . Because the general behavior of the Mrk 766 continuum is 
different from those observed from Sy 2, we expect that if these in- 
termediate velocity clouds are present in the NLR of Mrk 766, the 
fraction of the observed emission coming from these clouds should 
differs from the fraction in Seyfert 2. 

Another important point is that the flux in the IR region is 
higher for Mrk 766 than for Ark 564, which could indicate that the 
dust-to-gas ratio, d/g, is higher. Notice that d/g depends on the ratio 
between the dust re-radiatio n flux intensity and the bremsstrahlung 
in the optical - NIR range fcontini et al."2004'). As discussed in 
the following sections, observations constraints that ratio, and the 
requirement for high d/g becomes less strong. 

In summary, from the comparison of the SED of Mrk 766 and 
Ark 564, we conclude that both objects share a similar continuum 
energy distribution properties. It means that a safe starting point to 
model the line and continuum emission of the former is to assume a 
similar set of clouds successfully used in the modeling of the latter 
(Paper I) 



6 MODELING THE LINE SPECTRA 

By definition, the FWHM of the broad components of the per- 
mitted lines in NLSls are ^ 2000 kms^^ ( Osterbrock & Pogge 
Il985 h. Because the broad lines are narrower than those of classi- 
cal Seyfert 1 galaxies, it is often difficult to unravel the narrow and 
broad components of permitted lines. This usually leads to an in- 
correct determination of the flux contribution of the hydrogen lines 
emitted by the NLR, which in turn could induce large errors when 
comparing the observed and calculated line ratios. In order to cir- 
cumvent this problem, we chose to normalize the line intensities 
relative to that of [O III] A5007-I-4959 emission line flux as is shown 
in Table|2| 



The criteria adopted to model Mrk 766 follow the same line of 
reasoning adopted for Ark 564 (Paper I) and is summarized as fol- 
lows: (1) explain the emission-line spectrum and fit the continuum 
SED; (2) the velocities should be in agreement with the FWHM of 
the line profiles; (3) the ratio of calculated and observed lines rel- 
ative to H/3, R(H/3)mc, must be larger than unity, considering that 
the observed H/3 flux contains also the contribution of the BLR. 

A grid of models calculated for Ark 564 (Table 2 of Paper I) 
is rather complete, covering a large range of conditions suitable to 
describe the NLR of Seyfert galaxies. Moreover, recently, another 
specific grid of models was used to discuss the infrared continuum 
of Seyfert galaxies (Contini et al. 2004). All these results will be 
used as a starting point to explain the emission-line spectrum of 
Mrk 766. 

In the previous section, we saw that the continuum SED of 
Mrk 766 and Ark 564 are similar. It means that a starting point in 
explaining the emission-line spectrum of Mrk 766 should be the use 
of a multi-cloud model similar to that adopted for Ark 564. How- 
ever, as suggested by the FWHM of the observed line profiles in 
Mrk 766, the high Vs clouds found in Ark 564 (Paper I, models 
ml 3 and ml4 with Vs=1000 kms~^and Vs=1500 km s^^, respec- 
tively) are not expected to appreciably contribute to the observed 
spectrum. Thus, a first multi-cloud model is obtained for Mrk 766 
adopting the same single-clouds used in the Ark 564 but with a dif- 
ferent set of relative weights, which does not include models ml 3 
and ml4. 

In addition, a r ecent analysis of the optical - infrared contin- 
uum spectra of AGN. IContini et aljE004l) showed that low Vs (100 
kms^^) - high density (1000 cm~^) clouds are present in all types 
of Seyfert galaxies as well as clouds with intermediate velocity 
and density (Vs = 500 kms^^ and no = 300 s^'^). Since the ob- 
served continuum of Mrk 766 can not a priori exclude these type 
of clouds from the NLR modeling, a more consistent approach can 
be obtained by including them (models ml5 and ml6. Table |4}. 
The relative weights of the single clouds were changed in order 
to obtain a better fit to the observed emission-line spectrum. The 
physical properties of the single-cloud models, their H/3 absolute 
flux and the relative weights corresponding to the selected multi- 
cloud model (MC) are listed in the first rows of Table |4| (top). Ta- 
ble |4| also lists the contribution (in percent) of every single-cloud 
model to each emission-line. The (-I-) symbol near the line wave- 
lengths indicates that the doublet intensities are summed up. We 
have removed from Table |4| the contributions smaller than 3% in 
order to have a clear view of the results. The comparison between 
the observed and modeled emission-line ratios is made in Table|5| 
and|6l 

The multi-cloud model MC was calculated assuming that 
the radiation flux from the active centre is a power-law (F oc 
ly") of intensity Fh at 1 Rydberg (in units of photons cm~^ s~^ 
eV~^). The spectral indices auv = -1.5 (between 13.6 eV ^ 
E < 260 eV) and ax = -0.4 (for E> 260 eV) were obtained 
from the best matching model since the spectrum of Mrk 766 
is so rich in emission lines and continuum observations that it 
allows us to constrain the input continuum that produces them. 
This approach is preferred over the one that directly takes the ob- 
served far UV and X-ray data as input SED because the contin- 
uum that the NLR sees may differ from that emitted by active nu- 
cleus. In fact, othe r authors ( e.g. Sako et al. 2003; Page et al. 20(^ 
IVaughan & Fabia n 2003; Constantin & Shields 2003) have found 
spectral indices steeper than those used here for the ionizing radi- 
ation. In order to test this hypothesis, we have calculated single- 
cloud models assuming an ionizing spectrum with a power law 
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Table 4. Percent contribution of each model to the multi-cloud model MC 



d/g5 
wr 

[Nev] 3425 
[On] 3727+ 
[FeVII] 3758 
[Nelll] 3869+ 
[Sll] 4070 
[Olll] 4363 
Hen 4686 
[ArlV] 4740 
[Olll] 5007+ 
He: 5876 
[Fevil] 6086 
[Ol] 6300 
[Fex] 6375 
[Nil] 6548+ 
[Sll] 6717 
[Sll] 6731 
[Arm] 7135 
[On] 7320+ 
[Sin] 0.95 
[Sviii] 0.99 
[Sll] 1.028 
Hel 1.08 
[Six] 1.25 
[Six] 1.43 
[Sivi] 1.96 
Hel 2.06 
[Mgvill] 3.03 
[Siix] 3.93 



index Q(7v=-2.4 and ax = -1-7, which correspond to the indices 
reported by Constantin & Shield^ 12003). These new models were 
obtained with the same Vs and no as shown in Tableland were 
used to built the multi-cloud model labeled MCS. 

Although a reasonable fit to the data can be obtained with 
the cloud models listed in Table |7| we found that a high fraction 
of the line fluxes, corresponding to both high and low -ionization 
lines, are originated from high velocity clouds (Vs=1000 kms^^, 
model ml3 of Paper I). Such type of clouds are excluded because 
the measured FWHM of the line profiles in the nuclear spectrum 
of Mrk766 rules out the existence of gas motion with velocities 
above 500 km s^^ within the NLR. We therefore adopt model MC 
discussed above to obtain the SED of the continuum and to com- 
pare its prediction to the observations. Notice that the observed soft 
X-ray index is - 1 .77, while for the modeling we adopt -0.4. The dis- 
crepancy between both indices suggests that the NLR clouds see a 
harder, different continuum, than that we detect. This same conclu- 
sion had previously been reach ed by several authors jSinette et alj 
ll989lll993l:rKorista et alll997l) in studying the physical conditions 
of the BLR using photoionization models. IVIoreover, in our mod- 
els, the observed soft X-ray contiuum is due to bremsstrahlung 



orginating in the clouds that re-process the central ionizing radi- 
ation. 

Before concluding this section, it is important to comment on 
the values listed in Table|5| where the ratio between the calculated 
and observed line ratios relative to [O III], R([0 III])mc. and to H/9, 
R(H/3)mc, are listed. 

The calculated emission-line ratios for most lines except 
[IVIg VIII] and [S IX] are, within the error bars, good by a factor 
of 2. The latter two lines are underestimated by a larger factor. 

It is found that a large contribution to the bulk of the line 
fluxes comes from radiation-dominated clouds. This is the case, 
for example, of the NIR coronal lines. In them, the radiation flux 
from the active centre strongly contributes to the ionization of gas, 
which is heated to high temperatures (Tc > 1.5 10''' K) by the 
shock. Note that even in radiation-dominated clouds the shock is 
also accounted for. On clouds propagating outwards from the AC, 
it acts on the edge opposed to that facing the radiation flux. The 
coupled effect of the shock and the radiation flux determine the 
distribution of the temperature throughout the cloud. The temper- 
ature in the postshock region, immediately after the shock front is 
Te = 1.5 X lO'' (Vs/100)^ K decreasing downstream due to energy 
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Table 5. Comparison of observational data with model results. The indices MC and MCS identify the solution obtained with the two different continua used 
in the modeling. F or model MC, the input conti nuum is the one from the best matching model {ceuy = -1.5 and ax = -0.4). For MCS, the input continuum is 
the one derived bv lConstantin & Shield j<20O3ft . Se text for further details. 
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Table 6. Observed and predicted UV emission lines for the model MC. 
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losses by free-free, free-bound, line emission, and collision heating 
of dust. The stratification of the ions do vynstream follows the dis - 
tribution of the temperature described in Conti ni & Viegaj<199ll) . 
Our results suggest that the coronal lines in Mrk766 are formed 
mainly from the interaction of matter and radiation from the cen- 
tral source rather than shocks. A similar result is found from the 
analysis of optical and NIR coronal lines in a sample of Seyfert 1 
and 2 galaxies by Rodrfguez-Ardila et al. (2004, in preparation). 



They report, using high-spatial resolution spectroscopy, extended 
coronal emission within the few tens of parsecs from the active nu- 
cleus and points out that radiation from the central source are, in 
fac t, able to power them. Support for this scenario is also provided 
by jPrieto e t al. 2004) by means of VLT NIR observations of Circi- 
nus. Other evidence, however, supports the hypothesis that coronal 
lines are most likely associated to o utflowing winds( Erkens et alJ 
ll997tlRodriguez-Ardila et alj2002lj) . 
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Table 7. Percentage contribution of the individual clouds to the multi-cloud 
model MCS. calculated with steep spectral indices. 
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It is also possible to have a strong contribution from 
shock-dominated clouds, as is, for example, the ratio 
[O Il]3727-l-/[0 III]5007-I-. The high value of the ratio comes 
from the fact that 70% of the [Oil] 3727+ line flux comes from 
that type of clouds (models m8, m9, mil, and ml6). 

The optical [N Il]-l- lines are underestimated by the multi-cloud 
model MC by about a factor of 2. Since both R([Oin])i\/c and 
R(H/3)j\/c are less than unity, we conclude that the N/H abun- 
dance is about a factor of 2 higher than solar. 

On the other hand, all Fe lines are overestimated, indicating 
that Fe is depleted, probably trapped into small grains. Also, [Ar IV] 
4740-1- are overestimated, perhaps due to a too high Ar / H relative 
abundance. We notice that in order to improve the results for the 
Fe lines, both model MC and MCS require that Fe/ H be depleted 
by a factor of ~3. Furthermore, ratios of calculated to observed line 
ratios (relative to H/3), indicate an average contribution of the broad 
line region to the observed H/3 of about 40% (see columns 5 and 
7 of Table|5}. In Ark 564, that contribution was found to be lower 
(~30%). 

It is important to note that photoionization models for the 
NLR and the extended emission region of Mrk 766 were also pre- 
sented by GP96. They used two types of ionizing continual a 
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Figure 5. The observed and calc ulated continuum of Mrk 766 . The data 
come from NED (open circles) and lAlonso-Herrero et alJ |2223)- (solid cir- 
cles). The thin lines correspond to the results from the single clouds con- 
tributing to the best multi-cloud model MC, with dotted lines representing 
dust emission and solid lines bremsstrahlung emission from those clouds. 
The two vertical solid lines define the frequency range in which the X- 
ray radiation is absorbed. The spectrum obtained with the low velocity, 
high density single-cloud model ml6 is shown by the thin dashed line. The 
very thick solid line shows the continuum con'esponding to the multi-cloud 
model MC. 



power-law, with spectral indices -1 and -1.5 or a blackbody with 
effective temperatures rcff=10^ K and 2.2xlO^K. The compari- 
son between the calculated and observed emission line ratios of 
[O IIl]/H/3, [O Il]/H/3, [O I]/Ha, [N II]/Hq, [S Il]/Ha and [S IIl]/Ha 
showed that both the power-law and blackbody could explain the 
observed data. Clearly, the inclusion of emission lines located in 
other wavelengths intervals (NIR and UV), as was done here, gives 
more consistency to the modeling and allows a better understanding 
of the physical conditions in the NLR of Mrk 766. 

Finally, the observed and calculated UV lines relative to 
H/3 are compared in Table |6| The results show that there is a high 
contribution of model MC, i.e., of the NLR, to the UV lines, par- 
ticularly for N V, N IV] C IV and He II. 



7 THE CONTINUUM SED OF MRK766 

In order to consistently model the line and continuum spectra, it is 
necessary to cross check one another until a fine tune for both is 
found. Therefore, the results presented in the previous section were 
selected after this checking was performed. 

The calculated SED of Mrk 766 obtained from model MC is 
compared to the observations in Fig. El where the data take n from 
NED (open circles) and from Alons o-Herrero et alJ i2003h (solid 
circles) are plotted. Notice that the errors for the Alonso-Herrero et 
al.'s data are about 30% . The thin dotted lines and the thin solid 
lines represent the emission from the single clouds due to dust and 
bremsstrahlung, respectively. T he thin dashed line s correspond to 
the results from model ml6 (see lContini et alj2004l Sect. 3), which 
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Figure 6. The observed and calculated continuum of Mrk 766. The multi- 
cloud model was deiived assuming a relative weight for the low velocity, 
high density single-cloud (model ml6) a factor 10 higher than for MC (see 
text). Same notation as in Figurelsl 



is included in the multi-cloud model MC with other single-cloud 
models used to explain the Ark 564 spectra. The thick solid line 
represents the multi-cloud model with separated curves represent- 
ing re-radiation by dust and bremsstrahlung, while the very thick 
line represents the sum of these two curves. Finally, the dash-dotted 
line corresponds to the background old star population. 

Another SED, obtained by adopting a relative weight for 
model ml6 higher by a factor of 10 than the value used in model 
MC, is also shown in Fig|6| However, regarding the corresponding 
emission-line spectrum, this alternative multi-cloud model leads 
to very low line intensities relative to H/3 (R(H/3)a/c) and was 
dropped. 

The data in the radio range (long-dashed line) show that 
radiation at these frequencies is due to synchrotron created by 
the Fermi mechanism at the shock front. At long wavelengths, 
bremsstrahlung emission is self-absorbed, particularly for models 
with a high density downstream, which depends on no, Bo, and 
Vs. Actually, the optical thickness becomes larger than unity for 
T^ 10* K, ^ 10^^ Hz, for a slab of at least 1.3 pc downstream of 
high density clouds (e.g. model ml 2 ) because compression leads 
to densities higher than 2x 10 cm" <Osterbrock.l98&') . Thus, the 
contribution to bremsstrahlung in the radio range of relatively high 
density clouds must be negligible. 

Overall, the observed and calculated SED are in very good 
agreement, as can be seen in Figure|5| This confirms and strength- 
ens the validity of our approach. The largest discrepancy is located 
in the X-ray region, where the model overpredicts the observed 
emission. Here, the rapid and strong variability of the X-rays ra- 
diation, of up to a factor of 2 on a time scale of a few thousand 
seconds (Mason et al. 2003), clearly indicates that the bulk of this 
radiation is directly emitted by the central source and the BLR. 
However, it is expected that a small fraction of the X-ray emis- 
sion arises from high velocity clouds. Indeed, the bremsstrahlung 



calculated by models mil and ml2 overpredict the data. Notice 
however that a possible absorption of the X-rays by these kinds 
of clouds is not taken into account. In fact, absorption must oc- 
cur firstly because of the large geometrical thickness of the clouds 
(D=10 pc) with Vs=500 kms"^ and no=700 cm"^, and secondly, 
because of the high densities downstream due to compression (n/no 
10, depending on Vs, no, and Bo). Under these conditions, col- 
umn densities higher than 10^"^ cm~^ occur in the clouds producing 
a strong absorption of th e high energy radiation from 13.6 eV up to 
500 eV JZombeckl 1 990h . Clouds with lower values of D, Vsand no 
unlikely absorb X-rays. 



8 DISCUSSION AND CONCLUDING REMARKS 

As part of a wide program to understand the mechanisms at work in 
the NLR of NLS 1 galaxies, we have carried out a detailed study of 
the galaxy Mrk 766. For this purpose, the first spectroscopic obser- 
vations covering simultaneously the interval 0.8-4 /im is presented 
along with the measurement of the most conspicuous emission line 
fluxes detected in that region. These lines are used as additional 
constraints in order to unveil the physical conditions of the NLR 
gas. The single-cloud models adopted for modeling are those pre- 
viously used for the analysis of the continuum and emission-line 
spectra of Ark 564. However, the two high velocity cloud mod- 
els used for the later object were removed from the multi-cloud 
model that fits Mrk 766 because the X-rays flux is lower when com- 
pared to the continuum emission in other frequency ranges and the 
FWHM of the emission lines indicate lower velocities. On the other 
hand, two cloud models suggested by the fit of the NIR data of 
lAlonso-Herrero et ali (2003) sample were added. Indeed, after in- 
cluding these single-cloud models, a better fit of the line intensities 
relative to [O III] and H/3, as well as of the continuum, were be 
obtained. 

As already found in the analysis of other individual AGN, 
shock dominated models have higher relative weights because the 
absolute H/9 flux is lower. The shock velocities range between 100 
kms~^and 500 kms^^, the preshock densities between 100 cm~'^ 
and 1000 cm^'^ and the radiation fluxes from the active center from 
10^ to 5 10^^ photons cm^^ s^^ eV^^ at 1 Ryd with spectral in- 
dices of Q(7v=-1.5 and ax=-0.4. Higher power-law indices were 
dropped because they did not lead to satisfactory explanation of the 
observed spectra. 

The d/g ratios are between 10"^^ and 10"^* by number 
(4x10^^ to 4x10"* by mass for silicates) and reach the lowest 
value in clouds corresponding to model ml6 with d/g = 3 lO"^'^ 
(^ 10"^ by mass). 

We consider that R=Av/E(B-V) ~ 3.1 iCardelU et a^]|l989^ 
and Aa=1.086 t\ with rA=QA tt a-gr^ Ngr, where Q\ is the extinc- 
tion efficiency, ag^the grain radius, and Ngr is the column density 
of dust. Ngr = n X d/g x x, where n is the gas density downstream 
which increases with no and Vs, and x is the geometrical thickness 
of the dust slab. Adopting E(B-V)=0.5, 2, n=10*-10^ cm~^, 
d/g = lO"'^^ - 10~^*, X ~ 0.2 - 20 pc in good agreement with the 
geometrical thickness of the clouds given in Table 2. 

The important contribution to the NIR SED of Seyfert galax- 
ies of high density-low velocity clouds, corresponding to model 
ml6, was demonstrated by Contini et al. (2004), who also raised 
questions about their location in the galaxies and their dust-to-gas 
ratios. With the data already discussed, we are now able to roughly 
estimate the distance from the active center. To this purpose, we 
use the [O II] 7320-1- line, which is nicely fitted by the multi-cloud 
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model and shows the highest percentage contribution from model 
ml6 relative to the other forbidden lines. This type of cloud is also 
responsible for emission of the permitted lines He I 5876 A and He I 
1.083 /im. However, since the latter two lines include the contribu- 
tion from the BLR, any location of the clouds deduced from these 
lines may be mistaken. From Table |4] we see that the larger con- 
tribution to the [O II] 7320+ emission-line comes from the single- 
clouds mS, mlO and ml6. Thus, assuming that the single-clouds are 
co-spatial, their distance r from the central source can be estimated 
from the expression, 

Fyoiii] obsXd^ = {E,wr{i)xHl3{i)x{[OII]/HP),aUi))xr\4) 

where ;= 1,2,3, refers to the three types of single clouds, d is 
the distance to Mrk766, wr(i) and H/3(i) are given in Table|4] and 
the line intensity, relative to H/3, is equal to 4.64, 1.5, and 0.04, re- 
spectively, for m8, ml 1 e ml6. Assuming Ho=75 kms~^/Mpc, the 
distance of Mrk766 from Earth is d= 51.6 Mpc. Thus, we obtain 
an average distance of r = 160 pc. Since usually higher velocity 
clouds (m8 and ml 1) are closer to the center, r is probably a lower 
limit to the distance of the low-velocity clouds (ml6). The aver- 
age distance obtained for the clouds is also in very good agreement 
to the region size of 250 pc covered by the nuclear NIR spectrum 
employed in this work. This gives further consistency to our mod- 
eling an approach to describe the physical conditions of the NLR 
of Mrk766. 

The high no and the low dust-to-gas ratio of the clouds ml6 
may indicate that they are the residual of matter which was al- 
ready compressed by a previous blast wave propagating outwards. 
The velocity of the blast wave was high enough to destroy mos t of 
the du st grains by sputtering, i.e., VsJ^ 200 krns^^lContini et alj 
l2004h . Besides compression, another effect of the blast wave is 
fragmentation, which follows in the turbulent regime created by 
shocks. Indeed, in this model, because of the relative high density, 
the postshock cooling rate is high and the gas recombines at a dis- 
tance ^5 10^*^ cm from the shock front. Beyond this distance the 
gas temperature drops to < 500 K. Thus, the clouds are radiation- 
bound and the geometrical thickness (D) adopted for model ml6 
can be regarded as an upper limit. Note that the clouds correspond- 
ing to model ml6 could be geometrically thin, strengthening the 
hypothesis of fragmentation. 

In view of the important role of NLS 1 in the determination of 
metallicities in AGN, in particular through the nitrogen lines (e.g. 
[shemmer & Netzer 2002), the results of the multi-cloud model sug- 
gest that a N/H higher than solar by a factor ^ 2 would be in agree- 
ment with the [N Il]+/[0 III]+ line ratio (Table|5}. We also suggest 
that iron is depleted from the gaseous phase because of its emission 
line ratios are systematically overpredicted. 

Finally, let us conclude that the coupled effects of photoioniza- 
tion and shocks allowed us to explain most of the features observed 
in both the line and continuum spectra of the NLS 1 Mrk 766. 
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